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We study the charge state of the diffusing O2 molecule during silicon oxidation through hybrid 
functional calculations. We calculate charge transition levels of O2 in bulk Si02 and use theoretical 
band offsets to align these levels with respect to the Si band edges. To overcome the band-gap 
problem of semilocal density fuctionals, we employ hybrid functionals with both predefined and 
empirically adjusted mixing coefficients. We find that the charge transition level e^^~ in bulk Si02 
occurs at eV above the silicon conduction band edge, implying that the O2 molecule diffuses 

through the oxide in the neutral charge state. While interfacial effects concur to lower the charge 
transition level, our estimates suggest that the neutral charge state persists until silicon oxidation. 
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Determining the nature and origin of point defects at 
semiconductor-oxide interfaces is a major step in the de- 
velopment of electronic devices. - Since direct experimen- 
tal characterization of interfaces is often difficult, simula- 
tion methods based on density functional theory (DFT) 
are increasingly being used as an alternative tool for the 
investigation of atomic-scale properties. However, the 
study of defect levels within standard DFT schemes is 
still largely hindered by the well-known band-gap prob- 
lem. Indeed, calculated band gaps and band offsets are 
generally severely underestimated and the alignment of 
defect levels with respect to the relevant band edges is 
consequently affected. A correct alignment is a prereq- 
uisite for understanding issues such as oxide charging,^ 
defect passivation, Fermi-level pinning;^ stress-induced 
leakage current,'* charge state of defects, etc. 

Silicon oxidation is a fundamental process in which 
charging issues are expected to play a key vo\e^ Our un- 
derstanding of this process relies to a large extent on 
the model proposed by Deal and Grove, in which ox- 
ide growth proceeds through the diffusion of the oxygen 
molecule across the oxide followed by its incorporation at 
the silicon substrate.^ While supporting this general pic- 
ture, DFT calculations have been instrumental to achieve 
an atomic scale description of various aspects of silicon 
oxidation, such as the diffusion mechanism,—!^ the oxi- 
dation reaction)2ii£ the role of spiujii the layer-by-layer 
oxidation ) ^^1*'^ and the release of interstitial silicon..l4 
However, the charge state of the O2 molecule during the 
oxidation process has so far remained elusive. This is 
expected to determine the nature of the oxygen species 
occurring at the interfacei^ and consequently the oxi- 
dation reaction and the atomic structure at the Si-Si02 

Experimental attempts to determine the charge state 
of the oxidizing molecule rest on the study of the ox- 
ide growth kinetics in applied electric fields, but result 
in conflicting conclusions It should further be noted 
that the interpretation of such measurements is not triv- 



ial because of the complexity of the underlying atomic 
processes.^ Theoretical work based on gradient corrected 
DFT calculations concluded that the oxygen molecule 
diffuses through the oxide in a metastable neutral charge 
state, assuming the stable negatively charged state only 
in the vicinity of the substrate where electron tunneling 
can occurii^ This leads to a description of silicon oxi- 
dation in which neutral and negatively charged oxygen 
species play competing rolesjiS, This inference directly 
stems from the position of the charging level of the O2 
molecule relative to the silicon band edges. However, the 
theoretical determination of such an alignment is subject 
to the band-gap problem. 

In this work, we determine the charge state of the O2 
molecule during silicon oxidation by locating its charge 
transition level with respect to the relevant band edges at 
the Si-Si02 interface. To overcome the band-gap prob- 
lem, we use hybrid density functionals which give en- 
hanced band gaps compared to semilocal functionals. — 
Our results indicate that the neutral state of the oxygen 
molecule is thermodynamically most stable in bulk Si02 
for electron chemical potentials lying in the Si band gap. 
In the vicinity of the substrate, the image-charge inter- 
actions and the higher oxide density concur to lower the 
charge transition level. Nevertheless, our estimates of 
their effects suggest that the O2 molecules preserve their 
neutral state until the onset of oxidation. 

We considered the class of hybrid density function- 
als based on the generalized gradient approximation 
of Perdew, Burke, and Ernzerhof (PBE),*^ which are 
obtained by replacing a fraction a of PBE exchange 
with Hartree-Fock exchange. The functional defined by 
q;=0.25 is referred to as PBEO and has received some 
support from theoretical considerations.^'' In our calcu- 
lations, core-valence interactions were described through 
normconserving pseudopotentials and valence wave func- 
tions were expanded in a plane- wave basis set. The basis 
set is defined by an energy cutoff of 70 Ry. The pseu- 
dopotentials were generated at the PBE level and used 
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TABLE I: Band gaps of Si (Sf ) and Si02 and va- 

lence (AiSv) and conduction (AEc) band offsets at tlie Si-Si02 
interface, calculated at the PBE level, at the PBEO level, and 
through the use of the mixed scheme (Ref. 113). The experi- 
mental values are taken from Ref. |28|. 
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in all calculations. The electronic-structure calculations 
corresponding to large supercells were performed with a 
Brillouin-zone sampling restricted to the F point. For 
bulk silicon, we determined the band edges using a con- 
verged fc-point sampling. The integrable divergence of 
the Hartree-Fock exchange term was explicitly treated;^ 
Structural relaxations were carried out at the PBE level. 
We used the Quantum-ESPRESSC^^ and CPMI>^ codes. 

For the Si-Si02 interface, we used a superlattice struc- 
ture, in which a crystalline Si layer (9 monolayers) 
and an amorphous Si02 layer (~17 A) are periodically 
repeated.^"* The interface structure was designed to in- 
corporate a set of atomic-scale features inferred from ex- 
perimental data>^ The simulation cell in the transverse 
directions corresponds to a VS x -\/8 Si(lOO) interface 
unit. The model contains 131 Si atoms and 86 O atoms. 
For bulk amorphous Si02, we used a disordered 72-atom 
model structure obtained previouslyi^ 

First, we addressed the band alignment at the inter- 
face, as detailed in Ref. 127|. Band gaps calculated at the 
PBEO level are generally larger than the PBE values, but 
the agreement with experiment is not systematically im- 
proved (Table[l|. We determined band offsets by aligning 
bulk band extrema through a local reference level deter- 
mined in the interface calculation.^^ For crystalline Si, 
we took the average electrostatic potential as reference, 
but used the semicore O 2s levels for amorphous Si02 to 
overcome the difficulties associated with the structural 
disorder. As seen in Table |T1 the valence and conduction 
band offsets calculated at the PBE level underestimate 
the experimental results by about 2 and 1 eV, respec- 
tively. The use of PBEO, reduces these discrepancies by 
a factor of 2. To achieve a better description of the band 
alignment, we also considered a mixed scheme in which, 
for each interface component, we tuned the Hartree-Fock 
exchange fraction a to reproduce the experimental band 
gapi21 This resulted in a=0.11 for Si and q;=0.34 for 
Si02. The consistency of this scheme stems from the 
weak dependence of the interfacial dipoles on The 
band offsets obtained in the mixed scheme show excellent 
agreement with experiment (Table |T|. 

Second, we determined the relevant charge transition 
level between the neutral and the negatively charged 
state of the interstitial O2 in bulk Si02»^i^ The neutral 
O2 is repelled by the oxide network and is found at the 



FIG. 1: (Color online) Charge transition level £°' of the 
interstitial O2 molecule in bulk Si02 aligned with respect to 
the Si band edges, as calculated through the PBE, the PBEO, 
and the mixed scheme. Values in bold give the position of 
e''''" relative to the Si conduction band edge. Other values 
correspond to the Si band gap, the valence band offset, and 
the conduction band offset. Energies are given in eV. 

center of interstitial voids. ^ To account for the structural 
disorder in Si02, we considered 10 different interstitial 
configurations corresponding to representative local min- 
ima of the diffusion pathji The attaches to the oxide 
network by forming a bond to a Si atom.^^'^^ Our study 
comprises 12 of such atomic configurations. Formations 
energies of the oxygen molecule in its charge state q were 
determined as a function of electron chemical potential:^ 

E^i^,) = El, - - +qie. + ^, + AV), (1) 

where E^, , Ef^^^ , and are the total energies of the 
defect cell, of the bulk oxide, and of the isolated oxygen 
molecule, respectively. The electron chemical potential fi 
is referred to the top of the Si02 valence band £v, and AV 
is the correction needed to align the electrostatic poten- 
tial far from the defect to that in the unperturbed bulk 
(0.15 eV in our case). The total energies of the negatively 
charged species were corrected for the spurious interac- 
tion due to the periodic boundary conditionSfSi We set 
the charge transition level £^/~ in correspondence of the 
value of the electron chemical potential for which the av- 
erage formation energies of the neutral and the negatively 
charged oxygen species are equal.-?. Our PBE calcula- 
tions yielded at 4.2 eV from the valence band edge, 
in accord with previous results obtained with gradient- 
corrected functionalsjiii^ At the PBEO level, s^^~ lies at 
6.1 eV from the respective valence band maximum. In 
addition, we used a hybrid functional with Q!=0.34 in or- 
der to reproduce the Si02 band gap and obtained e^^^ 
at 6.6 eV from the corresponding valence band maxi- 
mum. The empirical adjustment of a is consistent with 
the mixed scheme for the band offsets. Furthermore, 
charge transition levels of the Si dangling bond obtained 
through such an adjustment of the band gap were found 
to agree with experiment within 0.06 eVr^ lending sup- 
port to this procedure. 

Third, we consistently aligned the calculated charge 
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FIG. 2: Schematic representation of formation energies of the 
neutral (bottom) and negatively charged (top) O2 molecule 
near the Si-Si02 interface as function of distance to the Si 
substrate. d-Si02 and Si02 represent the densified and the 
normal oxide, respectively. The wavy lines pictorially repre- 
sent the potential energy along the diffusion path. 

transition levels to the Si band edges through the theo- 
retical band offsets within each of the three considered 
schemes. The resulting alignment is shown in Fig. [TJ The 
three schemes give a similar picture, situating £^/~ be- 
tween 0.8 and 1.1 eV above the Si conduction band min- 
imum. Since the Fermi level at the Si-Si02 interface falls 
within the Si band gap, this result provides convincing 
evidence that the diffusing O2 molecule is thermodynam- 
ically most stable in its neutral charge state. 

We emphasize that this conclusion is naturally drawn 
provided the band offsets and the charge transition lev- 
els are calculated within the same theoretical scheme. 
For instance, the use of e°/~ calculated at the PBE level 
in combination with experimental offsets gives a charge 
transition level falling below the valence band maximum, 
leading to a qualitatively different description.^^ Such 
an alignment scheme assumes that the position of the 
charge transition level with respect to valence band max- 
imum is already well described in the semilocal DFT cal- 
culation. However, this is generally not the case. An 
improved description of the band gap results in signifi- 
cant shifts of both the valence and the conduction bands 
when referred to a reference potential external to the 
electron systemi^i^ At variance, charge transition levels 
of localized impurities remain largely unaffectedi^ This 
property suggests yet another alignment scheme not re- 
lying on the use of hybrid functionals. Assuming that 
the charge transition level calculated at the PBE level 
does not shift, we first determined valence band correc- 
tions relying on recent GW (—1.8 eV) and quasiparticle 
selfconsistent GW (—2.8 eV) calculationsi^ The respec- 
tive e*^/" levels are found at 0.5 eV and 1.5 eV above 
the silicon conduction band edge when aligned through 
experimental band offsets. This provides further support 
for the results in Fig. [TJ 
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FIG. 3: (Color online) Evolution of the charge transition level 
for the interstitial O2 molecule as a function of distance 
z to Si substrate. The charge transition level drops below the 
Si conduction band minimum at a distance zq ^\ K. We used 
esi =12 and esi02=4. The densified oxide (d-Si02, shaded) 
extends over ~ 10 A, beyond which the oxide (Si02) recovers 
its normal density. 



When the diffusing oxygen molecule approaches the Si 
substrate, two different interfacial effects acting on the 
neutral and the negatively charged O2 molecule need to 
be considered. X-ray reflectivity measurements indicate 
the occurrence of a densified interfacial oxide, showing 
a density of 2.4 g/cm^ and extending over a distance 
of about 10 The occurrence of a densified oxide af- 
fects the incorporation energy of the neutral O2 molecule 
which depends on the size of interstitial voidsi^ Conse- 
quently, the average formation energy of neutral O2 in 
the densified oxide increases by 0.6 eV.^''^^ This picture 
is confirmed by both selected calculations on our interface 
model and recent results in the literature f^^i^^ Figure [U 
schematically illustrates this effect showing the evolution 
of the formation energy of the neutral O2 as it approaches 
the interface. In the following, we express the correction 
with respect to the bulk formation energy as A£'''(z), 
where z is the distance to the substrate. 

The formation energy of the negatively charged oxy- 
gen molecule mainly results from the formation of a bond 
with a Si atom of the oxide network and is therefore not 
very sensitive to the modified structural properties of the 
near-interface oxide. However, as the negatively charged 
species approaches the silicon substrate, it undergoes sta- 
bilization because of the image-charge interaction arising 
from the dielectric discontinuity at the interface. From 
classical electrostatics, the stabilization energy reads 

Ai^-(.).--^ (2) 

where £si and £3102 a-re the static dielectric constants of 
Si and Si02, respectively. This effect is also illustrated 
in Fig. [1 

As a consequence of the two effects discussed above, 
the charge transition level also evolves as the O2 molecule 
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approaches the interface: 

= ^buik - ^E\z) + Ai?-(z), (3) 

where £^(,11^ corresponds to the charge transition level in 
bulk Si02. We remark that both effects concur to lower 
gO/-_ provide a quantitative estimate, we adopted for 
^buik the position determined within the mixed scheme, 
viz. at 1.1 eV above the silicon conduction band min- 
imum. The resulting evolution of e^^~{z) is shown in 
Fig. [3] We found that the charge transition level lies 
above the Si conduction band edge as long as the dis- 
tance of the O2 to the Si substrate exceeds ~1 A.— At 
this small distance the oxidation reaction is already un- 
der way^ suggesting that the negatively charged oxygen 
molecule plays only a minor role during silicon oxidation. 
The presence of charged oxygen species had previously 
been invoked to explain the growth kinetics in the thin- 



oxide regime However, alternative mechanisms based 
on a spatially varying diffusion rate do not require the 
charging of the diffusing O2 molecule.— 1^ In the present 
picture, the neutral O2 interstitial is the dominating ox- 
idizing agent. Other oxygen species can only occur upon 
the reaction of O2 with Si-Si bondsi^ 

In conclusion, we addressed the charge state of the O2 
molecule during silicon oxidation overcoming the band- 
gap problem through a scheme based on hybrid density 
functionals. The oxygen molecule is found to diffuse in 
its neutral charge state until the onset of oxidation. More 
generally, the defect level alignment scheme proposed in 
this work is expected to be very useful for addressing 
charging issues associated to defects at interfaces. 
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